Transport in Bilayer Graphene: Calculations within a self-consistent Born 

approximation 
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The transport properties of a bilayer graphene are studied theoretically within a self-consistent 
Born approximation. The electronic spectrum is composed of fc-linear dispersion in the low-energy 
region and fc-square dispersion as in an ordinary two-dimensional metal at high energy, leading to a 
crossover between different behaviors in the conductivity on changing the Fermi energy or disorder 
strengths. We find that the conductivity approaches 2e^/7r^/i per spin in the strong-disorder regime, 
independently of the short- or long-range disorder. 



I. INTRODUCTION 



Recently there was an experimental development in 
fabrication of atomically thin graphene, or single-layer 
graphite, that enables us to access its exotic electronic 
propertiesiSi^ The magnetotransport was measured and 
the integer quantum Hall effect was observed4s^ In the 
experiments, a multilayer that contains a few graphene 
sheets is also availableA^ The electronic structure of a 
bilayer graphene was studied theoretically and the spec- 
trum was found to be essentially different from that of a 
monolayer^ The purpose of this paper is to study trans- 
port properties of the bilayer graphene. 

A graphite monolayer has a /c-linear, massless Dirac- 
like spectrum and has long attracted theoretical interests 
as a "relativistic" problem in condensed matter physics, 
where k is the absolute value of the wave vector. The- 
oretical studies of transport in such an exotic electronic 
structure have been given by several authors, where the 
conductivity with or without magnetic field jSiS the Hall 
effect ^Siii quantum corrections to the conductivity,^^ 
and the dynamical transport^'^ are investigated. The re- 
sults show that the conductivity exhibits various singular 
behaviors in the vicinity of zero energyi2ii2iiiiiiiS4i& 

In bilayer graphene the energy dispersion includes both 
fc-linear and fc-square terms. In this paper, we calcu- 
late the diagonal conductivity in the absence of a mag- 
netic field within a self-consistent Born approximation 
(SCBA). A similar SCBA analysis has been applied for 
monolayer grapheneiSiiSiiSi We shall find that the coex- 
istence of fc-linear and fc-square dispersions leads to a 
crossover between transport properties similar to those 
of a monolayer graphene and to those of an ordinary two- 
dimensional metal as the Fermi energy is changed. The 
conductivity becomes nearly universal, 26^ /ir'^h per spin, 
in the case of large disorder. The analysis is made for two 
different kinds of scatterers, short range and long range, 
where the former represents on-site random energies dis- 
tributed on the carbon atoms and the latter a slowly 
varying random potential of the range much longer than 
the lattice constant but shorter than the typical electron 
wavelength. 

In Sec. II, the effective Hamiltonian and the resulting 
energy spectrum in a bilayer graphene are discussed, and 





FIG. f; (Left) Top view of the atomic structure in a bi- 
layer graphene. Solid and dashed lines represent the top 
[A and B sites) and bottom layers [A' and B' sites), re- 
spectively. (Right) Definition of the hopping parameters in 
a tight-binding model, 70 between nearest-neighbor sites in 
each layer, 71 between A and B' , and 73 between B and A' . 



model scatterers are introduced. The self-consistent Born 
approximation is briefly described in Sec. III. Explicit 
results are discussed for short-range scatterers in Sec. IV 
and for long-range scatterers in Sec. V. A discussion and 
a brief summary are given in Sec. VI. 



II. EFFECTIVE HAMILTONIAN AND ENERGY 
SPECTRUM 

The bilayer graphene is composed of a pair of hexag- 
onal networks of carbon atoms, which include A and B 
atoms on the top layer and A' and B' on the bottom, 
as schematically shown in Fig. ^ The two layers are ar- 
ranged in AB stacking, where A atoms are located above 
B' atoms, and B or A' atoms are above or below the cen- 
ter of hexagons in the other layers. The unit cell contains 
four atoms A, B, A' , and B' , and the Brillouin zone be- 
comes identical with that of the monolayer graphene. We 
model the system by a tight-binding Hamiltonian based 
on the Slonczewski-Weiss-McClure graphite modeliiiiS 
We include three parameters 70 , 71 , and 73 , where 70 rep- 
resents the intralayer coupling A B or A' ^ B' , and 
71 and 73 the interlayer coupling A B' and i? <-> A', 
respectively. The coupling parameters are estimated as 
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70 « 3.16 eV^ 71 « 0.39 eV^ and 73 « 0.315 eV.2i 

We can show that the low energy spectrum is given 
by the states around the K and K' points in the Bril- 
louin zone. Neighboring A and B' sites are coupled by 71 
to create the bonding and antibonding states away from 
the Fermi level, and the low energy states are given by 
the remaining A' and B sites »i The effective Hamiltonian 
reads, 
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where k± = k^ iky with k being the wave vector mea- 
sured from the K or K' points, and the effective mass 
TO* and the wave number /cq defined by 
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with the lattice constant a — 0.246 nm. The fc-linear 
term in the Hamiltonian ||2J) describes the direct hop- 
ping between A' and B sites, and the fc-square term the 
second-order process between A' and B via A-B' dimers. 
A typical energy where the fc-square and fc-linear terms 
become comparable can be defined by 
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The eigenenergy of Q becomes 
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where the upper sign corresponds to j = K and the lower 
to K\ s 
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fc^, and if 



arg(z) being the argument tp in z 
vectors corresponding to jnj are 




arg(fc+) with 
^. The eigen- 
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with 

Oku = ^Tgi-koe^'^ + ke-^^n 
Ok'u = arg(fcoe~''^-Hfce2"^). 



(8) 
(9) 



Figure 121 shows the energy dispersion © for j = K. 
In the high-energy region |e| > (l/4)£o, we have a single 
trigonally warped Fermi line, which becomes closer to 
a circle on going to the higher energy as the fc-square 
term dominates in the Hamiltonian. The equienergy line 




1 k^lko 



FIG. 2: Equienergy lines (top) and the three-dimensional plot 
(bottom) of the energy dispersion of the bilayer graphene 
around the K point. In the latter only the lower half (e < 0) 
is shown. 



becomes convex for e > lO.Seo- In the low-energy region 
\e\ < (l/4)eo, the Fermi line splits into four separate 
pockets for each of K and K', one center part and three 
satellite parts that are located trigonally. In the vicinity 
of zero energy \s\ <JC Eq, the dispersion becomes linear 
in k space with respect to the four Fermi points, where 
the center pockets can be approximated as a circle with 
radius k = je/eol^o and the three satellites as ellipses 
with the longer and shorter radii |e/eo|fco and |e/eo|A:o/3. 
While splitting of the four Fermi lines occurs only at very 
low energies, the trigonal warping extends to much higher 
energy as seen in Fig. |21 

The velocity operator = {dH/dkx)/fi has nonzero 
matrix elements between the states on the identical k 
points {j, k) written as 



{vjk)ss' = {j'ks'\vx\j'ks) = swjk + s'w*i^, (10) 



with 
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In the high-energy region |e| > (l/4)eo, the mean square 
velocity averaged on the contour e ~ |ejks| is given by 

(|(«,kWP>e = j^ + ^(l + (13) 

where the first and second terms come from the fc-square 
and fc-linear terms in the dispersion, respectively. In the 
vicinity of zero energy |e| <C eo, we have 
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The density of states per spin is defined by 

where J7 is the area of the system. In the high energy 
region |e| > (l/4)£o, this becomes 



Po(e) 
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with terms of the order of Oie^je^) neglected. In the 
vicinity of zero energy |e| ^ £oj o^i the other hand, 
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The density of states diverges logarithmically at |£|/£o — 
1 /4 due to the presence of saddle points in the dispersion. 

For the parameters mentioned above, m*/mo — 0.033 
with mo being the free-electron mass, A:o(27r/a)~^ = 
2.2 X 10"'^, and £o«3.9 meV. The electron concentration 
corresponding to (l/4)£o is = 1.7 x 10"'^° cm~^, that to 
£o is 1.0x10" cm-2, and that to IO.Seq is 1.1 xlO^^ cm"^ 
As the typical electron concentration in the present sys- 
tem is 10^^ cm~^ (£/£o ~ 10) or larger^ the trigonal 
warping is appreciable, but it is extremely hard to re- 
alize the situation where four Fermi lines are well split 
from each other. 

The Hamiltonian jSjl is formally equivalent to that of 
the monolayer system, where the /c-square term comes 
from a higher-order term in the k • p approximationi^ 
However, the parameter £o for the monolayer becomes 
the order of 70 and thus is larger than in the bilayer 
by the order of 1000, and the fc-square term gives only a 
small perturbation. In the bilayer, on the contrary, the fc- 
square term becomes dominant in the dispersion at small 
energy Sq. Note also that the Hamiltonian Q becomes 
invalid when the energy becomes as high as the anti- 
bonding states of the A-B' dimers. The deviation from 
(PI seems to become appreciable around£ ~ 71 /4 (~ 0.1 
eV corresponding ton^ ^ 2.7 x 10^^ cm~^)i The dis- 
persion of the bilayer graphene is much closer to that of 
three-dimensional graphite, where we see a similar trig- 
onal structure in the — ky plane with a fixed wave 
number along the stacking direction^L^ 



In terms of the eigenvector (t)^^" the amplitudes of the 
atomic orbitals at sites R"^ and are given by 

i^A'{^^') = ^0^^'^^exp[z(K,+k) -R^'], (18) 
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V'b(R^) = ^^i]'^exp[z(K,+k).R^], (19) 
V A* 



for j ^ K and AT', where ujk ~ ^jJk' = ^ with 

uj = exp(27rj/3), and N is the number of unit cells in the 
system. 

The dominant scatterers in the present system are not 
well known. In the following, we shall consider two kinds 
of model scatterers, short-range scatterers localized only 
on B sites or on A' sites and long-range scatterers, the po- 
tential of which spreads over a certain length scale larger 
than the lattice constant a. In each case we assume that 
the disorder strength is weaker than the dimer coupling 
71, so that we need to consider only the potential on 
A' and B sites. The effective Hamiltonian for the dis- 
order potential can be derived similarly as in monolayer 
graphene, if we identify the A' and B sites in the bilayer 
with A and B in the monolayer i2^i2^ 

The matrix elements for the short-ranged potential are 
written as 



ifk's'lUflj-ks) 
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(20) 



^ !|l-e^(K.+k-K,,-k').Rf (^^.,^/k'.y^^.^^k.^ (2^) 

for scatterers at R^ and Rf , respectively, where uf 
and uf are the integrated intensities of the potential, 
i.e., the coefficient of the 5 potential. For long range 
disorder, we assume that the potential range is much 
larger than the lattice constant a but smaller than the 
typical wavelength 27r/fc with fc being the wave number 
from the K or K' point. Then we can neglect the matrix 
elements for intervalley scatterings between K and K', 
while those for intravalley can be written as 

(/k's'|f/,b-ks) =<5,yge*(k-''')■^^'(c/,^"'''^')t(/,^■'^^ (22) 

where Ui is the integrated intensity of the potential. 

We have another possibility for the long-range disorder 
potential, where each of the scatterers is effective only 
in one layer. This is modeled by a long-range potential 
that has amplitude only over either A' or B sites. The 
situation then becomes almost equivalent to the short- 
range case, the only difference being that K and K' are 
decoupled, reducing the self-energy by a factor of 2. 



III. SELF-CONSISTENT BORN 
APPROXIMATION 

In the self-consistent Born approximation, the self- 
energy of the disorder-averaged Green's function {Ga,a') 
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is given by^ 



E 



(23) 



with a = (iks), where ( ) represents the average over 
the impurity configurations. In the present system the 
conductivity becomes isotropic in spite of the presence 
of strong trigonal warping. It can be calculated by the 
Kubo formula, 



a{e) 
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Re Tr [{v^G"'v.,G^) - (i-,G«z;,G«)] 



(24) 

where ^ {e - H + ii))-^ and = (e - W - ii))-^ 
are the retarded and the advanced Green's functions, re- 
spectively, with Ti being the Hamiltonian including the 
disorder potential. This can be rewritten as 



-Re Tr [v^{G^)i^^ {G^) - z;.(G«)^f ^(G«)] , 

(25) 

with v^^ = Vx{e + iQ,e — iQ) and v^^ = Vx{e + 10,6 + 
satisfying 



i^{e,e') = + {UG{e)i^G{e')U). 



(26) 



In the SCBA, should be calculated in the ladder ap- 
proximation. In the above and hereafter we omit the 
summation over the spin degeneracy, so the actual con- 
ductivity should be multiplied by a factor of 2. 

In the case of the present model scatterers, the self- 
energy and therefore the averaged Green's function be- 
come diagonal with respect to the wave number and the 
band index. Further, the self-energy is independent of 
the wave number and the band index, and thus is deter- 
mined by the energy alone. As a result we have 



a, a' 

Ga{e) = G{e,eo 



e - S(£) - Ea ' 
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where S(£) is the self-energy. 



IV. SHORT-RANGE SCATTERERS 

For short-range scatterers, we assume that they are 
equally distributed to A and B sites with density nf ~ 
ni/2 and the identical mean square amplitude 



{{ntf) = {{nf?) = 
by 



Then, the self-energy is given 
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1^ 



(29) 



By substituting the summation over j and k with the 
integration in energy e' = le^ksl: we can rewrite this as 



riiU 



de'po{e')J2G{e,se'). 



(30) 



The density of states is given by 
p{e) = — \-y"lmGc,{e + iO) = ^ImI](e + iO). 
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i 
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(31) 

For the conductivity, we can show that the vertex cor- 
rection vanishes in the short-range scatterers, or = v^, 
and we obtain 

jkss' 

" 2- Jo 

xRe [G(£ + iO, ss')G{e - iO, s'e') 

- G{e + iO,se')G{e + iO,s'e')]. (32) 

In the second equality we have replaced the summation 
over j and k with the integral in e' — |£jks|- 

We first look at the Boltzmann limit by taking S ^ 0, 
which should be valid in the weak-disorder case satisfying 
\e\ ^ |S]|. By assuming that E is infinitesimal in H3U|I . 
we can calculate S explicitly as 



-i-Weo = — , (33) 

2 Poo 2t 



where r is the lifetime and the dimcnsionless parameter 
W is defined by 



W 
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(34) 



For sufficiently high energy, h/2T sa {■k/2)£qW showing 
that W ^ 1 characterizes a typical disorder strength that 
smears out the fine low-energy structure due to the fc- 
linear term. This W is the same as that defined for a 
monolayer graphene in Ref. and as in Ref. 0. 

For the conductivity, we take only the terms with 
G^G^ and s = s' as a dominant contribution in l|32|l . 
and obtain a familiar form. 



cr(e) « e^poT(w^)e. 



(35) 



{\{vj\i)ss\^)e is the Fermi-surface average of 



Here 

the diagonal matrix element of Vx and use has been made 
of the approximation 



^Y.G".{e)Gt{s)^''-^p,. 



(36) 



With the use of (|13|) the conductivity for the higher 
energy |e| > Eq becomes 



1 f\e 



(37) 



The term proportional to |e| results from the fc-square 
term of the dispersion and is rewritten as Poo\£\t /m* 
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FIG. 3: Density of states per spin calculated in the SCBA. 
Plots become identical for the short- and long-range scatter- 
ers. 



as in the usual two-dimensional metal. The term inde- 
pendent oi \e\ is a correction due to the trigonal warping 
caused by the fc-linear term, which never vanishes even 
as £ — !■ oo. In the low-energy region |e| <^ Sq, we have 



(38) 



The energy-independent conductivity is essentially the 
same as that of monolayer graphene.^ 

In the case of large disorder W > 1, the fine structure 
in the density of states disappears completely (see be- 
low) and therefore the imaginary part of the self-energy 
becomes independent of energy and the real part becomes 
negligible. Thus, we always have 
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S](e + iO) w -iT, 



P(e) ~ Po 



(39) 
(40) 

(41) 



We can calculate the conductivity H32|l by substituting 
the high-energy expansions H13|l and (|16|l . and obtain 



a(£)«4-[5(e) + l] + 4-i;. 



with 
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e r\ e 
— I — arctan — , 
T e r' 



(42) 



(43) 



where terms of the order of 0{1/W^) are neglected. The 
first term of (|42|) comes from the fc-square term in the 
dispersion and the second term is a correction due to the 
/c-linear term. 
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FIG. 4: Calculated SCBA conductivity (solid) and Boltzmann 
conductivity (dotted) per spin in the case of short-range dis- 
order, with upper and lower panels showing smaller and larger 
ly's. For the Boltzmann conductivity in the lower panel we 
used the expression (l.'i?^ which is valid for e > eo. In the 
lower panel the approximate result given by Eq. 1421 is also 
shown. 



The conductivity goes to the Boltzmann limit 1)3 7|) for 
e > r, while for e < T it becomes 



e 



1 

W 



(44) 



It is interesting that the conductivity becomes universal, 
i.e., a —> 2e'^/Tr^h, and never vanishes in the limit of 
the large disorder. At a rough estimate, we can derive 
this expression by putting the uncertainty relation e ~ F 
in H37|l . In Sec. VI, this universal conductivity will be 
reconsidered in terms of Einstein's relation. 
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The self-consistent equation (|30ll can be solved easily 
by a numerical iteration. Figure 13 shows the density 
of states calculated for several disorder strengths. We 
notice that the logarithmic divergence present in po (e) at 
e/eo ~ 1/4 is smeared out very easily, and the structure 
in the vicinity of e = due to the fc-linear dispersion 
disappears around W ^ and is almost unrecognizable 
already for W = 0.4. 

The corresponding plot for the conductivity is shown 
in Fig. 21 along with the Boltzmann limit. In weak disor- 
der ^ 1 we observe a sharp dip at zero energy. This is 
the analog of monolayer graphenejS showing that fc-linear 
dispersion around zero energy remains intact at small 
disorder. Further discussion of the asymptotic value at 
£ = will be given in Sec. VI. The Boltzmann conductiv- 
ity drops to zero at |e| = (l/4)eo, as the velocity vanishes 
at the saddle points in the dispersion, but this singularity 
disappears due to the finite density of states. Apart from 
this difference the results are almost the same as those 
in the Boltzmann limit. The conductivities for large Ws 
are plotted in a different scale in the lower panel in Fig. 01 
compared with the high-energy expression of the Boltz- 
mann limit H37(l and the approximate expression (|42|l . 
The curves deviate from the Boltzmann limit in the re- 
gion e < F and never fall below 2e^/(7r^?i) even in strong 
disorder. The analytic expression becomes valid already 
for 1/2. 



LONG-RANGE SCATTERERS 



For the long-range disorder, we consider scatterers 
with density Ui and mean square amplitude ((wi)^) = u^. 
The two valleys K and K' are now decoupled and the 
conductivity is written as the summation of their indi- 
vidual contributions. The self-energy is given by 



E(e) 



2VL 



(45) 



which differs from if^ for short -range disorder in that 
the prefactor is larger by a factor of 2 and that the sum- 
mation is taken only in one valley. However those two el- 
ements cancel and the self-energy becomes identical with 
(|30|l in this notation. We can omit the valley index j 
completely. 

The conductivity without the vertex correction, de- 
noted by cr", is equivalent to (|32|) . The vertex part \2J6\ 
is given by 



(jks'|5^(£,e')|jks) = (jks>:r|jks) 
hU B{e,e') 



= ), (46) 



4m* 1 - n(£,e') 
where 11 and B are dimensionless quantities defined by 



n(£,£') 
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Gjks(£)G'jks'(£') 



(47) 
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FIG. 5: Calculated SCBA conductivity (solid) and Boltzmann 
conductivity (dotted) per spin in the case of long-range dis- 
order, with upper and the lower panel showing smaller and 
larger M^'s, respectively. For the Boltzmann conductivity in 
the lower panel we used the expression in H52|l which is valid 
for e > eo- In the lower panel the approximate result given 
by Eq. is also shown. 



xGjk,(£)G,k.'(£'), (48) 



k as' 



with Wj\^ — Wj]i{hko/4:m*) . The conductivity correc- 
tion including the valley degeneracy is then written as 



(49) 



TT^n 2W"^ V 1 - n^-^ 1 - n-R-R 

where B"^ = B{e+iO, e~iO) and B"^ = B{e+iO, e+iO) 
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In the Boltzmann limit, only the term including B^^ 
and n^'^ is relevant. A straightforward calculation gives 



n 



B 



B.A 



1 

2' 



Re(il;,ke-*«^"))e 



(50) 
(51) 



We can easily show that {Ke{wjke~^^^)) ^ = 1 to obtain 
6a/ {e^ /n'^h) = 1/W in the limit e > Eq. In the case 
£ ^ £oi on the other hand, (Re(i()jke^'^^''))e — and 
therefore 6a / {e^ /-k'^K) = 1/4W. The total conductivity 
combined with ct° in H37|) or (|38|l becomes 



a{e) 



1 

(£<£o)- 



(52) 



In the case of large disorder > 1, we can derive the 
analytic expression similarly to the short-range case. By 
using the self-energy (p??^ we can show that Bra = 1, 
Brr — 0, and Ura — 1/2, giving the vertex correction 
6a{e) ~ {e^ /Ti'^h)/W independent of energy. The total 
conductivity including cr" in H42|l is written as 



a[e) 



e 



S{e) + 1] 



For £ < r, this reduces to 

o2 



a 



e 



2 + A 

w 



(53) 



(54) 



The conductivity is given by a universal value as in the 
short-range case in the limit of large disorder. 

We show in Fig. the conductivity numerically com- 
puted for the several W^'s. The difference from the result 
in the short-range case is appreciable only in the clean 
limit W 1, because the vertex correction gives only a 
shift of the order of {e'^/h)/W. The conductivity again 
approaches 2e^/(7r^fi,) in the region £ < F in strong dis- 
order. The approximate result H54|) is valid for W > 0.8. 



VI. DISCUSSION AND CONCLUSION 

In Sees. IV and V we have seen that the conductivity 
in the case of strong disorder W > 1 deviates greatly 
from the Boltzmann conductivity and converges to the 
order of 2e'^ /ir'^h. This can be understood in terms of 
Einstein's relation 



2 7-1* 

e pD 



(55) 



with the density of states p and the diffusion constant D* . 
The diffusion constant is written as D* = {'Vx)t, where 
{v"^) is the average of the squared velocity over states at 
the Fermi energy and r is a relaxation time related to F 
through the uncertainty relation F = h/2T. 



If we neglect the fc-linear term in the dispersion, we 



have 



e 



and therefore a 



sB r/m* with 



ng — Poo|e|- This is nothing but the Boltzmann conduc- 
tivity for |£| ^ F. In the energy range |£| < F, however, 
we have (u^) ~ T/m* because we take an average over 
the states Sa ^ F. Thus, the diffusion constant becomes 
D* ~ h/2m* independent of energy using Fr ~ h/2. 
Upon using p = poo = m* /irh^, the conductivity be- 
comes a ^ /h. This conductivity is universal and in- 
dependent of the band parameters and the strength of 
scattering. It is also independent of short- or long-range 
disorder. 

The situation becomes different in the weak-disorder 
limit PF ^ 1, where a at zero energy drops from the 
Boltzmann conductivity almost by the factor W . This 
behavior is essentially equivalent to that in monolayer 
graphene, where the conductivity drops from the Boltz- 
mann to the universal value jir^h in the vicinity of the 
zero energy and the near-singular drop was ascribed to a 
reduction of the effective density of states contributing to 
the conductivity^ The universal value was shown to be 
unaffected by magnetic fields;^ and similar near-singular 
behavior was shown to be present in various transport 



quantities 



13.14.15 



A universal conductivity at zero en- 



ergy was also reported in the square tight-binding lattice 
model with one- half flux,l^- which has fc-linear dispersion 
as well. 

This zero-energy conductivity can be explicitly esti- 
mated from H32|) with the density of states H17|) and the 
square velocity H14(l . The integral turns out to be in- 
dependent of the imaginary part of the self-energy and 
returns a universal value 



-(0) - 5^, 



(56) 



which is six times as large as the conductivity in the 
monolayer. The extra factor comes from the product of 
the density of states and the square velocity, 3|£|/(7r^^), 
which is larger than in the monolayer by a factor of 6 due 
to the existence of the satellite elliptic Fermi pockets. 

Recently transport properties of bilayer graphene were 
studied experimentallyii& The resistivity exhibits a promi- 
nent peak at £ « and decreases rapidly with the in- 
crease of the energy or the electron concentration. This 
dependence is explained by Eq. H42I) or H53|l (or the lower 
panel of Fig.QlorO qualitatively quite well with the dis- 
order parameter 1 < IF < 2, in the region ris < 2 x 10^^ 
cm^^. In particular, the observed peak resistivity ~ 6.5 
kil corresponds well to the present result for VF ~ 2 in 
the case of long-range disorder (dominant usually) and 
IF '--^ 1 in the case of short-range disorder. However, the 
observed resistivity seems to decrease much faster than 
that given by Eq. H42|l or (|53|l with a constant W for 
larger electron concentrations ti^ > 2 x 10^^ cm~^. This 
might suggest that the effective range of the scattering 
potential can be comparable to the the electron wave- 
length at these electron concentrations. 

We have studied the quantum transport in bilayer 
graphene in zero magnetic field with the self-consistent 



8 



Born approximation. The coexistence of the fc-hnear and 
fc-square dispersion is observed as two different behaviors 
in the Boltzmann conductivity, hke that in the usual two- 
dimensional metal and in a monolayer graphene with a 
massless Dirac spectrum, which are energetically sepa- 
rated. The conductivity in the SCBA deviates from the 
Boltzmann limit in the case of strong disorder W > 1, 
and converges to ^ 2e'^ /{n'^h) (per spin). The conduc- 
tivity in the long-range scatterers exhibits a qualitatively 
similar behavior to that in the short-range case, while 
the vertex correction gives a positive shift of the order of 
l/W. 
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